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ABSTRACT 

The multiwavelength spectrum of BL Lacertae during its July 1997 outburst is 
analyzed in terms of different variations of the homogeneous leptonic jet model for 
the production of high-energy radiation from blazars. We find that a two-component 
7-ray spectrum, consisting of a synchrotron self-Compton and an external Compton 
component, is required in order to yield an acceptable fit to the broadband spectrum. 
Our analysis indicates that in BL Lac, unlike other BL Lac objects, the broad emission 
line region plays an important role for the high-energy emission. Several alternative 
blazar jet models are briefly discussed. In the appendix, we describe the formalism in 
which the process of Comptonization of reprocessed accretion disk photons is treated 
in the previously developed blazar jet simulation code which we use. 



Subject headings: Gamma Rays: observations — galaxies: active — BL Lacertae 
Objects: individual (BL Lacertae) 



1. Introduction 

A dramatic gamma-ray outburst of BL Lacertae (z=0.069) during July 1997 was recently 



reported by Bloom et al. 1997. BL Lac's gamma-ray activity was coincident with bright and 
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variable emission from optical to X-rays ( |Noble et al. 1997 ; Mattox et al. 1996); however, the 
light curves in each waveband are too sparsely sampled to determine accurately whether they 
are correlated. In addition, the gamma-ray spectrum was significantly harder than in the lower 
intensity state, which is commonly seen for flaring blazars (Mukherjee et al. 199(:; Mukherjee e\ 
al. 1997). Multiwavelength flares coincident with gamma-ray flares have also been documented 



for many other objects (eg., 3C 279, [Wehrle et al. 1998| and 1406-076, [Wagner et al. 19951) . The 
observations during the flare will be discussed further in Section 2. 

Historically, BL Lacertae was the prototypical 'BL Lac Object', i.e., has quasar-like continuum 
properties, but does usually not show emission or absorption lines (the equivalent width is < 5 
A). Though, Ha emission lines have been detected in BL Lac itself during a period of several 
weeks in 1995 (Vermeulen et al. 1995; Corbett et al. 1996| ). BL Lacertae is also a member of the 
blazar class of objects, i.e. a radio-loud AGN which has bright and variable broadband continuum 
emission as well as relatively high polarization. About 50 blazars have been detected in high 



energy gamma-rays by the EGRET instrument, and 14 are BL Lac objects (Mukherjee et al. 



1997) 



In this paper we will examine the multiwavelength spectral energy distribution (SED; or plot 
of vFy) of BL Lac considering several emission mechanisms as the possible source of gamma-rays: 
synchrotron self-Compton scattering (SSC; phisselini et al. 1985 ; Bloom &: Marscher 1996| ), 
external inverse Compton scattering with an accretion disk as the source of soft photons (ECD; 
Dermer et al. 199^ ), and an external inverse Compton scattering of soft photons reprocessed in 



broad line region clouds (EGG; [Sikora et al. 1994 ). These mechanisms and the techniques used 
will be discussed further in Section 3. The modeling results and their implications for the nuclear 
environment of BL Lac are presented and discussed in Section 4. A discussion of the application of 
alternative — leptonic as well as hadronic — models to the broadband spectrum on BL Lacertae 
follows in Section 5. We summarize in Section 6. In the appendix, we discuss our numerical 
treatment of Gompton scattering of reprocessed accretion disk radiation, which was not included 
in the original version of the jet simulation code used for our modeling procedure. 



2. Observations 

In addition to the optical (V and R band) and EGRET gamma-ray observations of Bloom et 
al. (1997), we have included additional optical {BVRI; |Webb et al. 1998 ), radio. X-ray, and low 



energy gamma-ray observations from the time of the outburst. The radio observations at 2.25 and 
8.3 GHz are from the Green Bank Interferometer monitoring program ( Waltman et al. 1996D , the 



4, 8 and 14 GHz measurements are from the UMRAO monitoring program ( lAller et al. 1985| ). 



The 2-10 keV ASGA X-ray spectrum ( Makino et al. 1997 ) shown in Figs. 1-3 (thick dashed 



curves) are from continuous measurements of July 18.6 - 19.62 and yielded a photon spectral 
index of 1.44 it 0.01. The ASGA spectrum agrees well with the XTE PGA spectrum of Madejski 
et al. (1997a) from the same epoch, which is not plotted in the figures. Madejski et al. (1999) 
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show that the X-ray spectrum is harder than for earher epochs, and that it is variable on time 
scales of approximately one day. The OSSE observations from 0.15 - 1.0 MeV have been provided 
by the OSSE Team ( Prove &: Johnson 1997 ). 



TeV upper limits from CAT ( |Fleury 1997| ) and HEGRA ( [Aharonian et al. 19991 ) from 



observations in 1997 August are also added (filled square and circle, respectively). Although these 
measurements are not quite simultaneous to the remainder of the observations, they might give a 
rough estimate of the actual TeV flux at the time of the July 1997 outburst because the optical 
flux during the TeV observations was similar to the flux observed during the EGRET flare. 

A reddening correction of A-^ = 0.7 derived from E{B — V) = 0.3 from Vermeulen et al. 
(1995) and the extinction law of Cardelli et al. (1989) were applied to the optical data presented 
in Bloom et al. (1997). Webb et al. (1998) have used a similar correction {E[B — V)= 0.36). The 
values of the optical flux plotted in Figs. 1-3 are averages during the week of the gamma-ray flare. 
The optical flux exhibited short-term variability on a timescale of several hours. The error bars 
shown in the figures correspond to the dispersion of the optical fluxes during the observing period. 



3. Outline of the model 

We modelled the broadband spectrum of BL Lac with the pair plasma jet simulation code 
described in detail in Bottcher et al. (1997, hereafter BMS). In this section, we only give a short 
outline of the model and the assumptions. 

Our model is based on various versions of leptonic jet models for blazar emission. A blob of 
ultrarelativistic pair plasma is moving outward from the central accretion disc along a pre-existing 
cylindrical jet structure, with relativistic speed /?rc and bulk Lorentz factor F. In the following, 
primed quantities refer to the rest frame comoving with the relativistic plasma blob. At the time 
of injection into the jet at height Zi above the accretion disc, the pair plasma is assumed to have 
an isotropic momentum distribution in the comoving frame. The distribution functions at this 
time are power-laws, n'{'y') = no 7'"*, with low- and high-energy cut-offs 71 < 7' < 72- The 
total density of electrons (and positrons) in the comoving frame is n'^. The blob is spherical in 
the comoving frame with radius R'^ which does not change along the jet. A randomly oriented 
magnetic field of strength B' ^ B'^p is present, where B'^^ is the magnetic field corresponding to 
equipartition with the energy density of pairs at the base of the jet. The jet is inclined at an angle 
Bobs with respect to the line of sight. The spherical geometry in the comoving rest frame implies 
that, in the stationary frame of the AGN, the blob appears Lorentz contracted with a length 
parallel to the jet axis of = R'^/T and a perpendicular extent R± = R'^. 

As the blob moves out, various radiation and cooling mechanisms are at work: synchrotron 
radiation, synchrotron self-Compton radiation (SSC; [Ghisselini et al. 1985 , Marscher Sz Gear 1985| , 



Maraschi et al. 1992| , [Bloom fc Marscher 199^ ), inverse-Compton scattering of external radiation 



from the accretion disc, either entering the jet directly (ECD for External Compton scattering of 
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direct Disc radiation; permer et al. 1992 , Dermer &: Schlickeiser 1992) or after being rescattered 



by surrounding (BLR) material (ECC for External Compton scattering of radiation from Clouds; 
[Sikora et al. 199^ , [Blandford fc Levinson 1995 , Dermer et al. 1997). The latter process was not 



included in the original version of the EMS simulation code. Therefore, its treatment is described 
in more detail in the appendix. 

The BLR is represented by a spherical layer of material with uniform density and total 
Thomson depth tt,blr-, extending between the radii Vin^BLR < ?" < i'out,BLR from the central black 
hole. The synchrotron mirror mechanism ( Ghisellini fc Madau 1996 , Bednarek 1998| ) is unlikely 



to be efficient in BL Lac objects ( [Bottcher &: Dermer 19981) , because it requires a high Thomson 



depth of the surrounding material {tt^blr ^ 0.1), for which no evidence exists in this class of 
objects. Therefore we do not consider it in our spectral modeling procedure for BL Lac. 

The central accretion disc is radiating a standard Shakura-Syunyaev (1973) disc spectrum 
around a black hole of 10^ Mg solar masses. The total disc luminosity is 10^^ L44 erg s~^. The full 
angle dependence of the disc radiation field, according to its radial temperature structure, is taken 
into account. Compton scattering of all radiation fields is calculated using the full Klein-Nishina 
cross section. The SSC mechanism is calculated to arbitrarily high scattering order. 

After the initial injection of ultrarelativistic pairs, all of the mechanisms mentioned above 
(except the synchrotron mirror) are taken into account in order to follow the evolution of the 
electron/positron distribution functions and the time-dependent broadband emission as the blob 
moves out along the jet. 77 absorption intrinsic to the source and the re-injection of pairs due 
to 77 pair production is taken into account completely self-consistently using the exact, analytic 
solution of Bottcher & Schlickeiser (1997) for the pair injection spectrum due to 77 absorption. 
However, with the parameters adopted to fit the spectrum of BL Lac, the blob is optically thin 
to 77 absorption for photons up to at least 100 GeV so that this process plays a minor role. 
Consequently, since we assume that the blob does not expand, the total electron (positron) density 
rig in the blob remains basically constant along the jet. 

The broadband spectrum investigated in this paper is a time average over the entire duration 
of the 7-ray flare. Thus, it corresponds to the time average of the radiation emitted by one or 
multiple plasma blobs rather than to a snapshot during the evolution of a single blob. To compare 
our model spectra to the observed broadband spectrum of BL Lacertae, we use the time-integrated 
broadband emission (fluence) of a single blob. An average flux is calculated by dividing the fluence 
by a representative time scale which corresponds to the time scale of repeated blob ejection events. 
A jet quasi-continuously filled with relativistic pair plasma corresponds to a repetition time of 
trep = 2-R^/(r c). We define a filling factor / of the jet given by / = f^"^^ /trep- In all our model fits 
we find that the repetition time scale (in the comoving frame) is generally longer than the cooling 
time scale for the highest-energy electrons, but comparable to the decay time scale of radiation 
produced by the inefficiently cooled part of the electron spectrum. Therefore, the variability in 
the high-energy portions of both the synchrotron and the inverse-Compton components (optical 
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— UV and > 100 MeV 7-rays, respectively) is dominated by light-travel time effects within the 
source. These effects have recently been studied in detail by Chiaberge & Ghisellini (1998). Strong 
variability on the light-travel time scale through the blob is expected in the high-energy parts of 
both radiation components, if the filling factor is substantially less than 1. The variability in these 
two energy bands (optical/UV and high-energy 7-rays) is expected to be correlated. If / ~ 1, 
any variability pattern is determined by variations of the electron injection parameters (i. e. the 
electron density, cutoff energies, spectral index etc.) and its modelling would require additional 
model assumptions about the underlying accretion flow and the pair injection mechanism. The 
low-frcqucncy portions of both components (radio — infrared and X-rays, respectively), which are 
predominantly produced by basically uncooled electrons, are expected to show relatively weak 
variability on longer time scales, also determined by long-term variations of the pair injection 
parameters. 

We point out that the evolution of individual blobs is virtually independent of the existence 
of multiple blobs along the jet if there are no significant fluctuations of the bulk Lorentz factor F 
within the flow. The radiation emitted by each individual blob is quasi-isotropic in the comoving 

frame so that the radiative interaction between different blobs may be neglected. The situation 
investigated here is geometrically very similar to shocks propagating along the jet because for 
large values of T the blobs are strongly Lorentz contracted along the jet axis. A major difference 
between a shock-in-jet model and our impulsive-acceleration model is that in the case of shocks 
particle acceleration will continue to take place at large distances from the central engine, where 
the external radiation field due to direct or reprocessed accretion disk radiation is negligible. 
Therefore, those models will be strongly SSC dominated. The spectral signatures will be basically 
indistinguishable from the SSC dominated cases of our model calculations, which we will discuss 
in the following section. 

Throughout this paper, we assume Hq = 65 km s^^ Mpc~^ and qq = 0.5. 



Apart from the observed broadband spectrum, there are several other observational results on 
BL Lac, which constrain the choice of parameters to be adopted in our modeling efforts. Recently, 
Denn et al. (1999) have observed superluminal motion in BL Lac of single components with up 
to P±^app. = (5-0 =b 0.2) = 7.7 lb 0.3. This indicates that the typical bulk Lorentz factor is 
expected to be in the range F ^ 8. This is much higher than the limit implied by the required 
transperency of the emitting region for high-energy 7-ray photons due to 77 absorption within the 
source, which yields 



4. Modeling results and discussion 




(1) 
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for a variability time scale of 1 day. Here D = (T [1 — Pr cos 9obs]) ^ is the Doppler factor. 

Occasionally, BL Lac shows strong emission lines. On May 21 and June 1, 1995, their width 
has been measured to be (v) ~ 4000 km s~^ (Vermeulen et al. 199£, Corbett et al. 1996| ). If we 
attribute this width to Keplerian motion, we find for the average distance of the BLR clouds from 
the central black hole: 



rBLR = 2.77-10-^M6pc. (2) 

We have used the XSTAR code (Kallman & Krolik 1998) in order to estimate the emission 
line luminosity resulting from a BLR consisting of Thomson thick (ry ~ 10) individual clouds, 
illuminated by disk blackbody radiation from an accretion disc around a black hole of Mg ~ 2, 
radiating at a luminosity of L^) ~ 10^^ erg s~^. We find that an angle-average Thomson depth 
tt,blr ~ 0.02 is sufficient to produce the observed line luminosity L^a ~ 10'^^ erg s^^. However, 
this value of the required tt,blr is strongly dependent on the dumpiness of the BLR and may 
well deviate by more than an order of magnitude from the value given above. 

The total luminosity of the accretion disc is poorly constrained since the spectral coverage is 
insufficient to clearly exclude the existence of a big blue bump from the accretion disc. We adopt 
a disc luminosity of = 10^^ erg s~^, implying that no big blue bump should be observable 
above the synchrotron spectrum. 

The shortest observed variability timescale during the 1997 flare was several hours, observed 
at optical wavelengths. The variability timescale at 7-ray energies was ~ 1 day, while it was 
slightly shorter than a day in X-rays. The optical variability yields a constraint on the extent of 
the emission region of 

R'b 5.4 • 10^^ D Ats cm (3) 

where Ats is the variability timescale in units of 5 hours. 

We have first attempted to fit the broadband spectrum of BL Lac with the basic models, 
where one of the radiation mechanisms mentioned above dominates the entire X-ray and 7-ray 
spectrum. 

Madejski et al. (1999) have presented an analytical estimate of the required parameters for 
an SSC fit to BL Lacertae, indicating that, if the entire 7-ray spectrum, peaking at esse ^ 2 • 10^, 
where e = hu / {m^c?) is the dimensionless photon energy, were produced by the SSC mechanism, 
unrealistic values of F ^ 100 and B' <^ 10~^ G would be needed. As mentioned in the previous 
section, the same arguments would apply to an SSC dominated shock-in-jet model. 

Fig. |l] illustrates our attempts to fit a strongly SSC-dominated model to the broadband 
spectrum of BL Lac. The optical flux and the hard X-ray and soft 7-ray spectrum are well 
reproduced by this model. It is also very well consistent with the concave spectral shape at 
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soft X-rays for which Madejski et al. (1999) have found weak evidence in ASCA observations of 
BL Lacertae in 1995 November, suggesting that the high-energy tail of the synchrotron component 
and the low-energy tail of the inverse-Compton component overlap in the 2-10 keV energy 
band. However, a pure SSC model is not able to reproduce the hard power-law above 100 MeV 
(EGRET). This is primarily because an SSC spectrum does not produce a strong spectral break 
in the MeV range which would be necessary to connect the X-ray / soft 7-ray spectrum to the 
EGRET spectrum. Our model 7-ray spectrum peaks at a values of esse significantly lower than 
2 • 10^, which is the only way to allow for more realistic values of F and B' than estimated above. 

In contrast to the SSC model, a pure external inverse-Compton model does produce 
the necessary spectral break in the MeV energy range and therefore reproduces the OSSE - 
EGRET 7-ray spectrum of BL lac better (see Fig. 2). However, the spectrum predicted by the 
ECD-dominated model is inconsistent with the ASCA X-ray spectrum. The jet filling factor 
corresponding to the simulation shown in Fig. ^ is / ~ 1 , which implies that any variability would 
be related to variations of the particle injection/acceleration mechanism. In order to have the 
ECD component strongly dominate the SSC component we need to assume a much larger mass of 
the central black hole (Mg = 100) than for our best fit presented in Fig. ^, where Mg = 2. 

The parameters required for an synchrotron + ECD model to reproduce the observed 
luminosities and peak frequencies in the broadband spectrum of BL Lacertae can be estimated 
using 

Lecd _ u'^cD _ 2Lji 
Lsy ^ u'^ ^ B'^zfcT^' ^ ' 

where Zi is the injection height of the plasma blob above the accretion disk, which is approximated 
by a point source located at the center of the AGN. Using 

B' ^ B,r^^^ ^0.9r-^G, (5) 
F escD 

where Bcr = 4.414 • 10^^ G, we find that eq. becomes independent of F, and 



Zi ~ J 1.6 • 10 pc, (6) 

Vcncr V ^ECD ^D^sy 

where ~ 10~^ is the dimensionless peak energy of the accretion disk spectrum. 

Since in the ECC process the accretion disc photon energies are boosted to much higher 
energies into the blob rest frame, because most of the soft photons enter the blob from the front, 
it is not possible to produce the OSSE - EGRET spectrum of BL Lac with a single ECC radiation 
component. However, it is appropriate to produce the observed 7 radiation in the EGRET 
energy range and, in particular, the spectral peak at eECC ^ 2 • 10^ with reasonable parameters. 
Approximating the external photon energy density due to reprocessed accretion disk photons as 
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4vr rsLR c 
we find 

Lecc ^ "^Lp Tt,BLR ( ^ECC 

or 



UeCC ^ A^^___2„ ' V'j 



(8) 



-2.3- 10^3 cm2, (9) 

tt,blr 

in reasonable agreement with our estimates on Tblr and tt^blr as given at the beginning of this 
section, considering the uncertainty of those estimates due to the unknown degree of dumpiness 
of the BLR. The co-moving magnetic field may be estimated as 



B' ^ Bcr^^eoT ^G.STG. (10) 

f-ECC 

The problems with fitting a one-component X-ray / 7-ray spectrum to the SED of BL Lac 
lead us to the conclusion that a second component is necessary to explain its peculiar high-energy 
spectrum. This is in agreement with the general tendency among different subclasses of blazars 
found by Ghisellini et al. (1998) that the sequence HBL LBL FSRQ may be related to an 
increasing dominance of external Compton scattering over synchrotron self-Compton scattering 
in the 7-ray regime. Since the properties of BL Lacertae are found to be intermediate between 
classical BL Lac objects and quasars, such a combination of SSC and ERG mechanisms appears to 
be a natural explanation for its peculiar high-energy spectrum. For the case of BL Lacertae this 
conclusion was first pointed out on the basis of analytical estimates by Madejski et al. (1997b, 
1999). Since the SSC mechanism yields a reasonable fit to the optical to soft 7-ray spectrum, we 
started out from this model and adjusted the parameters in a way that a second component shows 
up in the EGRET energy range. The fact that there is clear evidence for the existence of BLR 
clouds rather close to the central black hole, and the analytical estimates regarding this process 
suggest that the ECC mechanism is responsible for this additional high-energy component. 

Fig. m shows our best fit to the SED of BL Lac. The broadband spectrum from optical to 
high-energy 7-rays is reasonably well reproduced by our model. However, as all of our model 
calculations, it underpredicts the radio flux, because the low-energy cutoff in the electron energy 
distribution leads to a sharp cutoff in the synchrotron spectrum below ~ 10^^ Hz. Since there is 
no clear evidence for a correlation between radio and 7-ray flares in blazars ( Miicke et al. 19*9^ ), 



it seems plausible that most of the radio flux is produced in the outer jet regions, where the bulk 
kinetic energy of the jet is dissipated in interactions with external material. This is supported by 
the results of long-term radio, IR, optical, and X-ray monitoring of BL Lacertae as reported by 
Bregman et al. (1990). They have found that radio flux variations typically lack fluctuations at 
optical and IR frequencies by ~ 1 - 4 years. Their result that radio and IR flux variations are 
correlated with X-ray variations, but not with fluctuations of the optical flux, is consistent with 
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the IR and the X-ray flux being produced by strongly cooled electrons in the jet via synchrotron 
and SSC radiation, respectively. The radio, IR, and X-ray portions of the photon spectrum reflect 
long-term variations of the physical conditions at the acceleration site. Short-term variations of the 
injected electron spectrum are smoothed out by the slow cooling process. In contrast, the optical 
flux, produced by electrons cooling on much shorter time scales, reflects variations of the physical 
conditions on time scales comparable to the synchrotron cooling time scale of ultrarelativistic 
electrons. Using similar arguments, McHardy et al. (1999) have recently interpreted simultaneous 
IR and X-ray fluctuations of 3C 273 as evidence for the SSC process being responsible for the 
X-ray emission of that object. 

The fact that in Figs. || and ^ the total spectrum (solid line) is below the (unabsorbed) SSC 
spectrum (long dashed), indicates that here 77 absorption becomes important. As outlined above, 
77 pair production is included in our simulations self-consistently. 

The blob repetition time required by our fit implies a filling factor of / = 0.8, which means 
that the jet is filled essentially continuously. As explained in the previous section, this implies that 
due to light-travel time effects within the blob no strong variability on time scales shorter than 
the light-travel time through the blob are expected ( phiaberge fc: Ghisellini 1998 ), although the 



cooling timescale of ultrarelativistic electrons in the blob rest frame is only a few 100 s. Source 
variability is directly related to variations of the injection/acceleration mechanism. 

Due to the very time-consuming nature of our simulations, we cannot perform a detailed 
variational analysis of all parameters and give relative errors on their values. However, we point 
out that the relative strength of the different radiation components as well as their location in 
energy space constrain most of the parameters listed in the figure captions rather tightly. For 
example, an increase in the black hole mass leads to a) the BLR moving outward (Eq. [1]) and 
thus to a decrease of the ECC component, and b) the Schwarzschild radius and therefore the 
disc moving outward and becoming cooler, implying that the ECD component becomes stronger 
(because photons enter the blob under a more favorable angle), but emerges at different energies 
(combined effect of different disc temperature and different beaming pattern of disc radiation 
into the blob). Due to the colder disk, the ECC component will emerge at lower energies. Such 
changes also imply a change in the cooling time scale as compared to the required repetition time 
scale, determined by the total luminosity of a single blob which has consequences on the predicted 
variability pattern as mentioned in the previous section. This may indicate that the self-consistent 
system we simulate is rather sensitive to changes of the parameters, although there might still be 
ambiguities. 

Our results are sensitive to the assumed injection height Zi of the blob. According to the 
physical picture that the blob is produced by an explosive event in a hot corona above the accretion 
disk, we have generally chosen Zi ~ Rb, since there is no obvious mechanism how relativistic 
particles could be transported to significant distances above the disk without radiative losses. 



We find the appropriate average Thomson depth of the BLR, tt,blr = 0.025, to be 
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comparable to what we expect in other BL Lac objects (using the estimates given in Wandel 



1997). However, since the BLR seems to be located unusually close to the accretion disc 
(r = 5.5 • 10~^pc ~ 3 • Rs), it can still produce a significant flux in Ha line emission and 
contribute a significant soft photon density from rescattered accretion disc radiation. As pointed 
out in Bottcher &: Dermer (1998), a BLR of such low Thomson depth, located so close to the 
central black hole is extremely inefficient in terms of the synchrotron mirror mechanism, which 
provides additional justification for our neglect of this process. 

According to the results of Stecker &: de Jager (1997), 77 absorption by the intergalactic 
infrared background radiation is negligible for these photon energies at the redshift of BL Lac 
(z = 0.069), so the cutoff has to be intrinsic to the source. However, the CAT and HEGRA 
measurements have been performed in August 1997 and were therefore not simultaneous to the 
rest of the broadband spectrum. Since generally blazars show the strongest and most rapid 
variability at the highest photon energies, the > 300 GeV flux during the EGRET outburst could 
well have been higher than indicated by the upper limits in Figs. |l]-|^ 

Our model has generally problems with the hard optical spectrum. Although the correct 
optical flux is predicted, rapid cooling leads to a steep optical synchrotron spectrum oc z^"^/^ 
for a cooling particle population injected with a power-law 72(7) oc 7~^. The same problem 
has been encountered in other blazars, such as PKS 0528+134 (|Mukherjee et al. 1999| ) where, 
in some observing periods, the optical spectrum indicates an upturn compared to the overall 
synchrotron spectrum. This upturn is clearly inconsistent with the onset of the inverse-Compton 
component as well as with a "big blue bump" due to direct accretion disk radiation, which, in 
the case of PKS 0528+134, would imply an accretion disk luminosity of ~ 10^^ erg s^^. Up to 
now there is no generally accepted solution to this problem, and in most modeling efforts on the 
objects, it is simply ignored. A possible solution might be the effect of reacceleration, e. g. by 
hydromagnetic turbulences due to interactions of the jet material with the surrounding medium. 
Fermi acceleration by hydromagnetic turbulences will only affect the lowest-energy particles of 
the pair ensemble simulated here because for higher particle energies radiative losses strongly 
dominate reacceleration. Since the optical regime is close to the low-frequency turnover of the 
synchrotron component, reacceleration will thus have its dominant effect in this spectral range. 



5. Discussion of alternative models 

Until now, we have focused on the homogeneous, leptonic jet model characterized by 
instantaneous injection events of relativistic plasma blobs into the jet. We are now discussing 
a variety of alternative models which have also met with considerable success in explaining 
broadband spectra of blazars. For models invoking continuous injection of electrons into the jet, 
such as the accelerating-jet model discussed by Georganopoulos &: Marscher (1998), the final 
conclusions would be very similar to those found in the analysis of our model since our fits require 
jet filling factors close to unity. 
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5.1. Shock-in-jet models 

In the previous section we had already briefly mentioned that a shock propagating through 
the jet (e. g., Marscher fc Gear 1985, Kirk et al. 1998) would be strongly SSC dominated, and 



the broadband spectrum predicted by this model would be very similar to our SSC dominated 
case (see Fig. However, the spectrum could look drastically different if an external radiation 
component becomes efficient at large distances from the central engine. This external radiation 
component could most plausibly be provided by a cloud, reflecting part of the synchrotron 
radiation produced within the jet (the synchrotron mirror, Phiselhni fc Madau 1996| , [Bottcher fc| 



Dermer 199^). Bottcher & Dermer (1998) have shown that the energy density of the reflected 



synchrotron radiation in the comoving frame of the disturbance may be approximated as 



where Arc and tt^c are the radial extent and the Thomson depth of the reflecting cloud of cold 
material, and z is the distance of the disturbance from the central engine. In this model, R'^ 



may be identified with the extent of the disturbance along the jet (denoted x in Marscher & 



Gear 1985). Approximating further u' TT^dil) u'^, where rr,d is the Thomson depth of the 



disturbance along the jet axis, and (7) T ^Rsy/^sy, we find 



^Rsy _ . ^ ^ -p <^Rsy (. 2VR'g 

— 4 TT.d TT.r, i 1 



L., -4r^,.rT,er^^(l-^). (12) 



Eq. 12 indicates that the synchrotron mirror becomes efficient, if z ^ ^ ^'b- following, 
we will require this condition to be fulfilled. The variability time scale of a high-energy 7-ray flare 
in this scenario is given by At^^r = max {Arc/(r'^c) , -R^/(rc)}, which yields a constraint on Arc 
in addition to the constraint on R'^ as given in Eq. ^ Using these constraints, we find 

^ 1 LRSy ^sy o in-12 

T^T,dTT.c^-— ~ 8 • lU , [16) 

4 Lsy eRSy 

which is appears to be a reasonable number. Assuming that the hard X-ray / soft 7-ray spectrum 
of BL Lacertae is dominated by the SSC process — as we expect in the shock-in-jet model — , 
and that the SSC component peaks at 2 ^ esse ^ 200, we derive a bulk Lorentz factor of 
10 ^ r ^ 100. Then, the condition z » cT'^ At^ar yields ^ 1 pc. Since the disturbance is 
emitting synchrotron and SSC radiation basically independent of the external (reflected) radiation 
field, we would expect an extended high state of 

Atsj/,5SC ~ > Atvar,EGRET (14) 

at mm, infrared, optical, UV, X-rays, and soft 7-rays prior to the EGRET flare. This seems 
to contradict the observed quasi-simultaneous flaring at optical frequencies and in the EGRET 
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energy range as reported by Bloom et al. (1997), unless the synchrotron mirror operates at very 
low efficiency due to z ~ 2r R'^. 

5.2. Leptonic cascade models 

Blandford & Levinson (1995) have proposed a model invoking a pair cascade developing along 
the jet as synchrotron and inverse-Compton radiation produced by the accretion disk and within 
the jet is rescattered into the jet by circumnuclear material and initiates a pair cascade via 77 pair 
production with the inverse-Compton emission produced by relativistic electrons injected into the 
jet at the injection radius. In its original form, this model predicts single power-law high-energy 
spectra, but it is well conceivable that the soft and hard 7-ray spectra are composed of different 
radiation components, as in our model. One problem with this pair cascade model is that it 
predicts that the high-energy emission in the EGRET regime lags variations at lower frequencies, 
and the high-energy variability should be characterized by longer variability time scales, which 
is in contradiction to most 7-ray flares observed from blazars so far. This is because the 7-ray 
photosphere, beyond which the jet becomes optically thin to 77 pair production for photons of a 
given 7-ray energy, increases with increasing photon energy. Also, the observed quasi-simTiltancity 
of the optical and 7-ray flares of BL Lac with probably roughly equal time scales might be hard 
to reproduce with this model. 

A different model, also based on a leptonic pair cascade, has been suggested by Marcowith et 
al. (1995). In this model, the jet consists of a mildly relativistic, cylindrical outer electron-proton 
jet which powers an ultrarelativistic, inner pair jet via hydromagnetic interaction at the boundary 
between these two regions. A pair cascade is initiated by 77 absorption of inverse-Compton 
scattered direct accretion disk photons with themselves and the radiation field of the accretion 
disk. The model predicts a strong spectral break (a^ « 2 ax) at the photon energy where 77 
absorption becomes important. As in the Blandford & Levinson (1995) model, this scenario 
predicts a 7-ray photosphere increasing with 7-ray energy. The Marcowith et al. (1995) model 
has been very successful in explaining the high-energy spectrum of 3C 273, where many other 
models have problems with the large spectral break around ~ 10 MeV. However, in its published 
form, the model does not provide detailed predictions about the synchrotron component, and it 
predicts a smooth single power-law at high energies, inconsistent with the MeV - GeV spectrum 
of BL Lacertae. An additional soft photon source for inverse-Compton scattering could possibly 
solve the latter problem, but these issues would need to be addressed in detail by the authors of 
this model, before a final conclusion about its applicability to BL Lacertae can be drawn. 
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5.3. Hadronic models 



In this subsection, we discuss models in which most of the relativistic particle energy in 
the jet is primarily contained in protons. Their energy can either be transferred to secondary 
pairs through pj initiated pair cascades or via collisionless acceleration mechanisms to primary 
electrons. 

Mannheim et al. (1996) have presented a model fit of the proton-initiated cascade (PIC) model 
of Mannheim (1993) to a non-simultaneous broadband spectrum of BL Lacertae. This model seems 
to have the potential to provide a reasonable fit also to our simultaneous broadband spectrum. 
The hardening of the 7-ray spectrum at several GeV may be interpreted as the blue-shifted bump 
resulting from the decay of first-generation vr'' mesons produced in pj interactions. Flares at 7-ray 
energies are expected to be simultaneous to synchrotron flares, in agreement with the observations 
of BL Lacertae. For model parameters derived in the context of this model see Mannheim et al. 
(1996). 

The second class of initially proton-dominated jet models is based on plasmoid deceleration 
and collisionless transfer of internal energy of swept-up protons to non-thermal electrons ( |Pohl &: 
Schlickeiser 1999| ). The radiation signatures of these processes have been investigated in detail by 
Dermer & Chiang (1998) and Chiang & Dermer (1999), and their application to blazars has been 
discussed in Dermer (1999) and Pohl &; Schlickeiser (1999). For moderate (F ~ 10) bulk Lorentz 
factors of the plasmoid SSC is so far assumed to be the primary 7-ray production mechanism. 
The simultaneity between the synchrotron and the 7-ray flare observed in BL Lacertae suggests 
that the radiative cooling time scale of non-thermal electrons in the blast wave is shorter than the 
dynamical time scale of the system, i. e. the electrons are in the fast-cooling regime. In this case, 
and if SSC scattering occurs in the Thomson regime, the 7-ray luminosity may be approximated 
as 



Lssc 1 fp-2\ Ce , . 

(15) 



Lsy 3 VP - 1/ es 

where p is the spectral index of the nonthermal distribution into which electrons are accelerated, 
and ee and are the energies in nonthermal electrons and in the magnetic field as a fraction of 
their respective equipartition values. For a standard value of p ~ 2.5, this yields eg ~ 30 e^, which 
implies that the magnetic field strength has to be substantially below the equipartition value. This 
has also been found for 7-ray bursts by Chiang &: Dermer (1999), where too strong a magnetic field 
would lead to too rapid electron cooling, producing too soft X-ray spectra, inconsistent with the 
hard low-energy slopes observed by BATSE. A pure SSC spectrum, however, would not produce 
the peculiar spectral shape of the high-energy spectrum observed in BL Lac. A solution to this 
problem could be a blue-shifted vr'^ bump ( Pohl fc Schlickeiser 1999| ) or an additional contribution 
from Comptonization of an external radiation field. 



- 14 - 



6. Summary and conclusions 

We have analyzed the multiwavelength data for BL Lacertae during its July 1997 outburst 
in terms of several popular models for the broadband emission of blazars. We found that 
a reasonable fit is only possible with a combined SSC/ECC model. The parameters of the 
central accretion disc and its environment which yield an acceptable fit are Mbh = 2 • 10^ Mq, 
Lfl = 10^^ erg s~^, tj-^blr = 0.025. The minimum variability time scale consistent with our model 

is ^tuar^min ^ 1*2 h. 

The need for an external inverse-Compton component in the GeV photon energy range, 
in addition to the SSC component, confirms that BL Lac is an atypical BL Lac object. For 
example, the simultaneous broadband spectra of the X-ray selected BL Lac objects like Mrk 421 
(pVIastichiadis &: Kirk 199^ ) and Mrk 501 ( plan et al 1998| , |Petry et al 1999| ) can very well be 



fitted with a pure SSC model, while the emergence of an external inverse-Compton component 
seems to be more typical of quasars, where strong line emission indicates the existence of relatively 
dense broad-line regions. In fact, most simultaneous broadband spectra of quasars during 



outbursts can well be fit with an ECC or ECD model (e. g., 3C273: Dermer et al. 1997 ), or a 



combined SSC/ECD/ECC model (e. g., PKS 0528-^134: [Bottcher fc Collmar 19981 [Mukherjee ei 



al. 19991 ; 3C279: [Hartman et al. 1999^ . This trend has been more firmly established by a detailed 



modeling analysis of 51 7-ray loud blazars by Ghisellini et al. (1998). 

Since BL Lac belongs to the subclass of radio selected BL Lac objects whose properties are 
generally intermediate between quasars and X-ray selected BL Lacs, our results might strengthen 
the hypothesis that it is the density and structure of matter surrounding the central engine of 
an AGN which determines its appearance as either quasar or BL Lac object. In particular, this 
implies that the differences between these classes of objects are not predominantly due to viewing 
angle effects. 

Our conclusions are certainly biased in the sense that we have concentrated on leptonic jet 
models of impulsive particle injection, which is the only type of models which we have supported 
by detailed numerical simulations. We have done so because of the considerable success with which 
this model has met in explaining a wide variety of blazar spectra in many different 7-ray intensity 
states. We have discussed several alternative models, such as shock-in-jet models, leptonic cascade 
models and hadronic jet models. Several of these possibilities could not be ruled out. Detailed 
numerical calculations of the application of those models to our simultaneous broadband spectrum 
of BL Lacertae are beyond the scope of this paper, but we strongly encourage the authors of those 
models to do similar analyses so that a definitive conclusion in favor of a single model can be 
reached. 
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A. Reprocessed accretion disc radiation 

For the effects of Compton scattering of accretion-disc radiation rcscattcrcd into the jet 
trajectory by surrounding clouds, the central accretion disc may be treated as a point source, 
radiating isotropically, since in the process of rescattering in the BLR, virtually all spatial 
information about the point of origin of photons in the accretion disc gets lost. Furthermore, 
the accretion disc spectrum is approximated by a single-temperature blackbody spectrum. The 
blackbody temperature Qd = kTii/{mf,(?) is determined in a way that the resulting spectrum 
peaks at the same photon energy as the total photon spectrum of the Shakura-Sunyaev (1973) 
disc. This yields the photon production rate of the disc, 

where K is determined by the normalization to the total disc luminosity Lu, 

K = sr-^ = 1.88 . 10^9 ^ s-^sr-^ (A2) 

Throughout the appendix, the asterisk denotes quantities in the accretion disc rest frame, whereas 
quantities without asterisk denote quantities in the blob rest frame, e = hvl{me(?) is the 
dimensionless photon energy, and $7^ = (/x|j, (ffj^) is the solid angle of emission of an accretion disc 
photon with respect to the jet axis. 

Following the treatment of rescattering in surrounding material as outlined in Bottcher & 
Dermer (1995), but for a time-independent situation, we find for the differential photon number of 
rescattered accretion disc photons at height z above the disc, in the stationary frame, 

n\ie* n*-z) - ' f dr fA3) 

n^^ie ,SZ ^g^2c^eve^_l J '^\\z^ - - rz f,*j,y ^^^^ 
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where r is the distance of the scattering location in the BLR from the center of the accretion disc 
(the black hole), and x = ^r'^ ^ — 2rzfi*j^. follows from simple geometrical considerations, 
ne(r) is the electron density in the BLR at the location r which is determined by r and ^*^. 
This photon density is transformed to the blob rest frame using the Lorentz invariance of 
np/,(e,f])/e2 = n;^(e*,f]*)/e*2and 



e* = er(l + /3r/u), (A4) 

where ji and ji* are the angle cosine under which the photon enters the blob in the comiving 
and the stationary frame, respectively. The photon spectrum (A3) is used in the procedure 
described in BMS to calculate the electron energy-loss rate and the emitted photon spectrum due 
to inverse-Compton scattering. 

If the blob is located well within the inner boundary of the BLR, z <^ rin,BLR, the energy- loss 
rate can be well approximated assuming that all rescattered photons enter the blob from the front, 
i. e. with /i = — 1. This is the head-on approximation. Furthermore, in this case we approximate 
the photon spectrum by a delta-function in photon energy. 



e3 



,^r(l+/3r/.) 



where 



B2 S{e - (£)) (A6) 



and 

" 15 Vr(i + /?r/^) 

Combining these two approximations yields 

S, head— on 



^,.l!f,^V. (AS, 



where (e) = 2.7T@d- The integral I involves the angle integrations over the Klein-Nishina cross 
section, to which an analytical solution is given in BMS. In our simulations, this approximation is 
used as long as z < 0.2 • rin^BLR- In this range, its deviation from the exact electron energy-loss 
rate is negligible for all electron energies. 

If 
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given in the blob rest frame by 




r(i+^rM) 



-Ir, 




oo 



r\z'^ — x"^ — r z 



ne(r)x 



(All) 



where e = 6^/(7^ [1 — Px]) and x = M/^s is the cosine of the collision angle between the scattering 
electron and the photon before scattering. Here, we have neglected the very weak dependence of 
the scattered photon spectrum on the azimuthal angles involved. In our simulations, we use Eq. 
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10^ io'' 10'^ 10'^ 10'' 10'^ 10'' 10'^ 10'^ 

v[Hz] 

Fig. 1. — Attempt to fit the SED of BL Lac witli an SSC-dominated model. Parameters: 71 = 3000, 
-y2 = 10^ s = 2.6, B' = 1.2 G, r = 15, Bobs = 2^ n'e = 10^ cm-3, i?^ = 10^5 cm, Zi = 10'^ pc, 
/ = 0.6 
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10^ io'' 10'^ 10'^ 10'' 10'^ 10'' 10'^ 10'^ 

v[Hz] 

Fig. 2. — Attempt to fit the SED of BL Lac with an ECD-dominated model. Parameters: 71 = 30, 
^2 = 10^, s = 2.6, B' = 0.8 G, r = 20, Oobs = 2", < = 400 cm'^, R's = 2- lO^^ ^m, Zi = 7-10-3 pc, 
Me = 100, / = 1 
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v[Hz] 

Fig. 3. — Fit to the SED of BL Lac with a combined SSC + ECC model. Parameters: 71 = 500, 
72 = 3.5 • 10'', s = 2.5, 5' = 9 G, r = 15, 9obs = 3°, < = 3 • 10^ cm-^, R'^ = 1.2 • lO^^ gm, 
Zi = 5- IQ-^ pc, rin,BLR = 5 • 10"-^ pc, rout,BLR = 6 • 10"^ pc, TT,BLR = 0.025, Me = 2, / = 0.8 



